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The possible existence of four-neutron resonances close to the physical energy region is explored. 
Faddeev-Yakubovsky equations have been solved in configuration space using realistic nucleon- 
nucleon interaction models. Complex Scaling and Analytical Continuation in the Coupling constant 
methods were used to follow the resonance pole trajectories, which emerge out of artificially bound 
tetraneutron states. The final pole positions for four-neutron states lie in the third energy quadrant 
with negative real energy parts and should thus not be physically observable. 
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1—5 1 I. INTRODUCTION 

oo : 

The existence of pure neutron systems could have far-reaching implication in nuclear physics 0. However the 
question of multineutron existence is far from being cut-and-dried both theoretically and in experiment, being fed 
' by a long series of controversial predictions and measurements. Recently much attention have been paid to the 
Cn| ■ possible existence of bound tetraneutron (four- neutron system). This interest has been triggered by the experimental 
observation of few events in the 14 Be break-up reaction 2] . On the other hand such a prospect rises serious objections 
from the point of view of nuclear interaction theory. It has been shown by several groups 0, 0, IE that realistic 
nuclear Hamiltonians exclude the existence of bound 3 n, 4 n and even larger neutron clusters. In fact, the most 
favorable mechanism to construct tetraneutron would be by putting together two virtual (almost bound) dineutron 
pairs. However, in order to force the binding of virtual dineutrons one has to have very strong neutron-neutron 
interaction in P and/or higher partial waves, which is not compatible with our present comprehension of the nuclear 
interaction. 

' Nevertheless the possible existence of resonant states in pure neutron systems having observable effects in nuclear 
reactions, could not be eliminated. Such a scenario is evoked in a recent analysis of 8 He(d, e Li)4n reaction: some 
excess of low energy 6 Li nuclei has been observed, which can not be explained by phase space analysis involving both 
four free neutrons and two non-correlated dineutron pairs in the final state ■ Furthermore authors of in their very 
recent study Q agreed that previously observed signals could result from the existence of near-threshold four-neutron 
resonance, without involving bound tetraneutron. The aim of this study is to clarify whether or not the existence of 
resonant tctraneutrons can be tolerated by modern nuclear interaction models and thus if these experimental claims 
can be supported in a theoretical ground. This work is a natural extension of our preceding work |(|, in which wc 
have demonstrated that realistic nuclear Hamiltonians exclude the existence of physically observable three-neutron 
resonances. 

No proper ab-initio calculations of the resonant tetraneutron with realistic nn forces are known to the authors. 
Some conclusions were drawn in favor of its existence based on calculations of tetraneutron bound in an external 
well ||; furthermore it was suggested that these resonances could have rather large widths. The only rigorous 
study of tetraneutron resonances was accomplished in ref. |lOl | using the simplistic MT I-III nn interaction, which 
contains only 5-waves. Unfortunately no observable resonances have been found there and only the existence of some 
broad subthreshold structures (S-matrix poles with negative real energy parts) was pointed out. The same authors 
remarked however that the positions of these subthreshold states strongly depend on the details of the nn interaction 
used. Realistic nucleon-nucleon (NN) models contain indeed interactions in higher partial waves and are therefore 
better suited to accomodate tetraneutron and push its resonant states out of the subthreshold region. 
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II. THEORETICAL BACKGROUND 



Although many of the nuclear excited states are resonances, they are seldom considered in theoretical nuclear 
structure calculations due to the huge technical difficulties of solving the continuum states in many-body problems. 
These states are often treated as being bound, but such a procedure is justified only for very narrow resonances and 
is not appropriate in our case. Resonant tetraneutron, if existing at all, will probably have a rather large width. The 
problem we are dealing with represents therefore a double challenge: first it is a four-particle problem and second, 
being a continuum state, it has an exponentially diverging wave function. We will present in what follows the equations 
allowing to solve the four-particle problem in a mathematically rigorous way and describe the methods that make 
possible the treatment of resonant states. 




FIG. 1: Faddeev-Yakubovsky components K and H . Asymptotically, as z — » oo, components K describe 3+1 particle channels, 
whereas components H contain asymptotic states of 2+2 channels. 

In order to solve the four-body problem we decompose the wave function into a sum of 18 Faddeev-Yakubovsky 
(FY) components, see Fig^ an d rewrite the Schrodinger equation as a set of coupled FY equations [TJ]. If all four 
particles are identical, only two of the 18 FY components are independent, which we denote by K and H. These 
components are furthermore related by two intcgrodiffcrcntial equations: 

(E-Ho-V)K = V(P+ +P-)[(1 + Q)K + H] 

(E-H -V)H = VP[(1 + Q)K + H] (1) 
where P + , P~ , P and Q are particle permutation operators: 

P + = (P-)- = P23P12; Q = eP 34 ; P = P13P24 = P24P13, (2) 

and £ is a Pauli factor related to the exchange of two identical particles: e = — 1 for two identical fermions. Using 
these notations, the four-body wave functions is given by: 

* = [1 + (1 + P+ + P~)Q] (1 + P+ + P-)K + (1 + P+ + P")(l + P)H (3) 

Each FY component F = (K, H) is considered as a function of a proper Jacobi coordinate set (x, y, z) , defined 
respectively by 

xk = r 2 -ri _ _ 

/— - _, %H = ^2 — T\ 

VK - V^ 3-11 ^) VH = ? A -r 3 (4) 

z K = /f(r 4 -*+*+*) & = V2(^-^) 

The angular, spin and isospin dependence of these components is expanded using tripolar harmonics y a (x, y, z), i.c: 

(xyz\F) = J2 Fa{x ' y ' z) y a (x,y,z). (5) 

The quantities F a (x, y, z) are called regularized FY amplitudes, and the label a holds for the set of 10 intermediate 
quantum numbers describing a (J n ,T = 2,T 2 = 2) state. When describing tetraneutron, the isospin dependence of 
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the FY amplitudes can be omitted, and the set of quantum numbers a reduces to 8 elements. We use the j — j scheme 
for the intermediate coupling of FY amplitudes, defined as 



K = < (lx(si82)<T m ) jx (l y 8 3 )j y (lzSi)j,\ (6) 



(lx(siS2)<r x ) jx (i»(s3S4)<7 1) ) , . lz\ (7) 



H = 

where is the spin of the individual particle and the total angular momentum of the four-particle system. 
Each of the N c = Nk + Nh amplitudes in the expansion (JjjJ is further conditioned by the antisymmetry properties 
(— ) tT *+ i *+ 1 = e for K and (— Y x+lx+1 = (— ) a y +l v +1 = e for H. FY components K and H are regular at the origin, 
and it can be shown that for bound state problem, they decrease exponentially outside the interaction domain. In 
this case one can impose these functions to vanish on the borders of some constrained box: 

F a (x — Xmax, V — ^max, 2! = Zmax) — (8) 

Equations 111181) are enough to solve the bound state problem. 

Resonance wave functions are however divergent and cannot be described by the boundary conditions ijSJ. In 
order to solve the resonance problem, we make use of two different methods, successfully applied in [t| to treat the 
three-neutron system. The implementation of these techniques in the four-body FY equations is analogous to the 
three-body case. Therefore we only briefly discuss them here and the interested reader can refer to for technical 
aspects. 

The method of Analytical Continuation in the Coupling Constant (ACCC), proposed by Kukulin et al. [l^, is 
based on the fact that a resonant state arises from a bound one when the interaction between the particles is made 
less attractive. The corresponding eigenenergy is considered as an analytical function of a coupling constant A, which 
determines the strength of the attractive part of the potential. Therefore, one can try to analytically continue the 
energy of the bound state as a function of the strength A to the complex plane and obtain this way the width and 
the position of the resonance. It can be shown moreover that close to the threshold, where bound state turns into the 
resonance, the momenta k = \JE — Eq is proportional to 

A — An for virtual state 



VA — Aq for resonant state 



(9) 



where Ao is a critical value of the coupling constant and Eq = E(Xq) is the threshold energy. If multipartiele system 
does not possess bound states in its subsystems, as it is a case for multineutron, then Eq = E(Xq) = 0. 

It turns out that using an analytical continuation of k(x) in terms of a simple polynomial expansion converges 
slowly and that Pade expansion of order [N,M] 

un, m( ^_ a\x + a 2 x 1 + ... + a N x" 

( >~ l + bl x + b 2 x2 + ... + b M xM [W) 

is more appropriate. 

It is quite simple to put ACCC method in practice. One should artificially bind tetraneutron by adding some 
attractive interaction to the system's Hamiltonian H = Ho + W a tt ■ Then, the critical value of the coupling constant 
(Ao) is determined and several eigenenergies values E (A*) are calculated for Aj > Ao; they are used to fix the Pade 
expansion (|10f> coefficients. However, to make this extrapolation efficient, one should provide rather accurate binding 
energies values E (A,) and an especially precise Ao as an input. While only few lowest order terms in Pade expansion 
are enough to determine the positions of narrow nearthreshold resonances, the description of deep resonances requires 
several terms and very accurate input of E (Aj). The determination of high order Pade expansion terms requires at 
least five digit accuracy in the binding energies. 

The other method we use, namely Complex Scaling (CS) [Tflj . can be applied to calculate resonance positions 
directly. This method makes use of the similarity transform 

S = e l6r &, (11) 

applied to the Hamiltonian of the system, i.e: 

(SHS- 1 ) (S* res ) = E res (5* res ) . (12) 
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Such transformation does not affect the eigenvalue (E res ) spectra. However, if the scaling angle is large enough 

9 > I |argi? res | - the modified resonance eigenfunctions (^S^ res ^j become square intcgrablc. Evidently, the CS 

method can be applied to FY equation. By this transformation all the radial variables r = (x,y,z) in cq.Q arc 
replaced by re 10 = (xe ,ye l9 , ze ). The problem becomes analogous to a bound state one with complex variable and 

transformed FY amplitudes {^F^j , which unlike resonance eigenfunctions ^> res arc in Hilbcrt space. 

CS transformation requires the analytical continuation of the potential V(x) into the complex plane V{xe lS ). This 
turns out to be a weak point of this method when applied to nuclear systems, since as discussed in nuclear 
potentials have mischievous analytical properties: they become strongly oscillating and even divergent already for 
relatively small transformation angles 8 > 30°. This fact limits the applicability of CS method to narrow resonances, 
with Im(—E res ) < 2Re(E res ) values. 

The numerical solution of FY equations is performed by expanding F a (xyz) on a basis of three-dimensional piecewise 
Hermite polynomials and projecting equation Q with boundary conditions (|SJ| onto tripolar harmonics. In this way, 
the integro-differential FY equations are converted into a linear algebra problem: 

AX = E res BX 

with A and B being large square matrices, whereas E res and X are respectively the eigenvalue and the eigenvectors 
to be determined. The reader interested in a detailed discussion on the formalism and the numerical methods used 
should refer to @- 



III. RESULTS AND DISCUSSION 



The results presented below have been obtained by using the charge-symetry break ing Re id 93 potential to describe 
the nn interaction. This choice is dictated by purely practical reasons: as discussed in [SJLlJIi Reid 93 model has better 
analytical properties to perform the complex scaling operation than his coordinate-space modern concurrents. 
We would like to remark however that other realistic NN interaction - namely AV14, AV18 and Nijm II - provide 
very similar results for two- and three-neutron systems and exhibit also a similar behavior for the - artificially bound 
- tetraneutron. These facts let us believe that no qualitative changes in the four-neutron resonance can emerge from 
the properties of a particular model. All the calculations presented in what follows use the value ^— = 41.44 MeV-fm 2 
as an input for the neutron mass. 

In a similar way as in the study of three- neutron system we introduce an additional attractive four-nucleon 
(4N) force to analyze the tetraneutron resonance trajectories in a systematic way. We have chosen the form: 

V 4n = ~W pe~k , (13) 

where W and po are respectively the strength and range parameters of the potential, and the hyperadius p = 
\J x 2 + y 2 + z 2 is an invariant quantity with respect to the permutation operators 0. Such kind of force is easy to 
implement in FY equations In our previous work [t| devoted to the three-neutron system, we have used a 3n force 
having the standard Yukawa form. However we have found the functional form l|13(l more appropriate for studying 
artificially bound tetraneutron. This form does not diverge as p — > and thus avoids a rapid shrinking of the bound 
structures generated. 

As it has been already remarked in an extremely strong additional interaction is required to force tetra- 

neutron binding. As a consequence, the thus generated bound system is a very compact object making unlikely its 
physical existence. On the other hand, resonances are extended structures. In order to ease the transition from bound 
to resonant tetraneutron, we have fixed a rather large value for the range parameter po m (|13[l and taken po = 2.5 
fm, a value considerably larger than the one we could expect for a realistic 4N interaction. 

Our strategy to study 4n-resonances is to vary the strength of the potential W and trace the resonance energy- 
trajectory E res (W). The final resonance positions, which correspond to realistic nuclear interaction, are eventually 
reached at E res (W = 0). 

When applying ACCC method, the parameter A given in is identified to the 4N force strength A = W . We 
determine several auxiliary values of Ei(Xi) in the bound tetraneutron region A.; > Ao. These values are later used as 
an input to determine the Padc expansion coefficients (at, bj) of ea. Hl(J|) . Few calculations are performed in low energy 
region as well to determine the critical value Wo = Ao for which tetraneutron is bound with zero energy E res (Wo) = 0. 
Once (fli, bj) and Ao are known, we use equation (|10fl to analytically continue E res (W) curve to the resonance region 

w < w . 



5 



When applying CS method, we perform a series of direct resonance calculations for several decreasing values 
of W < Wq, until the calculations become unstable due both to the large size of the resonance widths and to the 
necessity of using ever increasing scaling parameter 6. This method is used in the near threshold region - W is Wq and 
Im(— E res ) <<Re(E res ) - where it gives very good results to improve the accuracy of Aq used in ACCC calculations. 
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TABLE I: Critical strengths Wo (MeV-fm^ 1 ) of the phenomenological force 11311 required to bind four neutrons in different 
states. The range parameter po of this force was fixed to 2.5 fm. W' are the strength values at which tetraneutron becomes 
subthreshold, i.e. Ke(E res ) = 0. In the last row of the table are given the physical resonances positions (W = 0). 

In Table [i] are summarized (second row) the critical strengths Wq required to bind tetraneutron in states with 
different J T quantum numbers. Even though we have taken 4N force with a rather long range, the critical strengths 
values Wq are still considerable. Other noticeable feature is that these critical values are almost equal for all negative 
parity tetraneutron states J n = — , 1 _ and 2~ we have considered. The reason for such a degeneracy is that 
tetraneutron binding energies are insensitive to nn interaction in P and higher partial waves. Actually, their values 
remain unchanged up to three-four digits if these l nn > 1 interaction terms are switched off. Tensor coupling is present 
only in l nn > 1 partial waves and has very small impact on these states. As a consequence the total spin (S=l) and 
angular momenta (L=l) are separately conserved. 

A similar situation is observed for positive parity states J T = + and 1 + , which are also almost degenerate. These 
states are dominated by the FY amplitude of type K with l x — l y — l z = intermediate quantum numbers, i.e. 
they are almost pure L=0 states. On the other hand, they differ by their total spin (S=0 for J v = + and S=l for 
J w = 1 + ). Unlike expected, the J™ = + state has a structure dominated by a "nn+n+n" type configuration and not 
by a "nn+nn" one, i.e. containing only a single 1 Sq dineutron pair and not two. For J 71 " = 2 + , Wq is considerably 
larger. This state must have a total spin S=2 to be realized with L=0, i.e. all neutron spins pointing in the same 
direction and thus no any 1 Sq dineutron pairs with antiparallel neutron spins is present. The corresponding large Wq 
value can therefore be understood as the price to pay for breaking the remaining dineutron pair. J n = + and 1 + 
states remain also unchanged if nn P-wave interaction is switched off. A very strong enhancement of these waves is 
required, as much as creating a dineutron resonance!, in order to see their effect in the binding energies. The only 
state sensible to nn P-waves is the 2 + . 

In figures |21 and |3 arc displayed the tetraneutron resonance trajectories for the same negative and positive parity 
states we have considered in Tabled I n both figures, CS results are indicated using empty symbols (square, cercle and 
triangle) which correspond to different values of the 4NF strength parameter W. ACCC trajectories for different states 
are depicted by solid, dashed and dotted lines. They have overimposcd star-like x,*,+ symbols which correspond to 
the same W values than those used in CS calculations. In order to compare quantitatively the agreement between 
both methods, the numerical values of some resonance positions are also given in Table ITT1 This agreement is rather 
nice for narrow resonances. For wider resonances, small discrepancies appear, which are due to the drawbacks present 
in CS method described above. Sizeable differences appear for resonances with Ke(E res ) < —Im(E res ), which is the 
limit of applicability of CS transformation. 

In this region, ACCC results are still rather well converged (better than 5%) with respect to the Pade expansion. 
However this convergence is getting worst when one departs further and further from the bound state region. The 
accuracy of Pade expansion is 20% near the subthreshold region, where resonance trajectory moves into the third 
energy quadrant. The accuracy of the physical resonance positions, when the additional interaction is fully removed 
(W = 0), is at worst of 50%. 

These limitations in ACCC accuracy are due to the increasing size of the Pade expansion argument x = y/X — Ac 
which forces to take into account higher order terms. The precise determination of high order Pade coefficients fails 
due to the severe accuracy criteria it imposes to the input. As an example, we have illustrated in figure |U the Pade 
expansion convergence for J w = 0~. We can see that the shapes for [N,M] = [3,3] and [4,4] order Pade expansion 
curves are already very close to each other. However the separation between the energies corresponding to the same 
W values still exist and it increases when one departs from the bound state region. 

The uncertainty in determining the final resonance positions is also manifested in figurc[5| In this figure, we compare 
the resonance trajectories for J 11 = 2~ tetraneutron state obtained with 4NF of ea. (|13fl having different values of 
the range parameter pq. For p = 2.5 fm we ended with a value E res = —1.4 — 9.7i MeV, while for pq = 2 fm with 
E res = -2.5 - 12.5* MeV. 
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FIG. 2: Negative parity tetraneutron resonance trajectories parametrized by the strength W of the phenomenological 4NF. 
ACCC results are denoted by lines with overimposed x, *, + symbols. They correspond to W values by steps of 4 MeV-fm" 1 
starting from 38 MeV-fm" 1 for J n — CP and 2" states and from 36 MeV-fm -1 for 1". CS results are represented by circles, 
squares and triangles. 




FIG. 3: The same as in figure |21 for positive parity states. W values are reduced from 20 MeV-fm 1 with step of 2 MeV-fm 
for J n = + , from 21 MeV-fm" 1 in step of 2 MeV-fm" 1 for 1+ and from 38 MeV-fm -1 in step of 4 MeV-fm" 1 for 2 + . 

As it has been discussed above, tetraneutron negative parity states on one hand and J" = + and 1 + ones on 
the other hand, are almost degenerate in energy. This degeneracy is also reflected in the corresponding resonance 
trajectories, which superimpose close to the threshold. Notice however that the small difference - not exceeding several 
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TABLE II: Comparison of CS and ACCC method results. Resonance positions for tetraneutron states obtained by adding 
phenomenological 4n force with strength W (in MeV-fm - ) and range po=2.5 fm are compared. 
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FIG. 4: Convergence of ACCC method with respect to the order [M,N] of Pade expansion for J" = 0~ tetraneutron state. 
ACCC curves are followed by star-like points indicating the resonance positions for W values decrasing from 36 MeV-fm -1 by 
step of 4 MeV-fm - . CS results are presented by full cercles and corresponding to W values from 36 to 15 MeV.fm. 



keV - in the binding energies, results into an increasing separation of these curves. This demonstrates the necessity 
of performing very accurate inputs for the Pade extrapolation and the difficulty of obtaining accurate predictions for 
broad resonances in the ACCC method. 

Regardless the convergence problems mentioned above, our results indicate that the final resonance positions will 
always stay in the third energy quadrant for all tetraneutron states. An accurate determination of the physical 
resonance position is not possible with the methods used in the present work. Nevertheless, for all calculations we 
have performed, they were situated in the third energy quadrant (Re(E)< 0,Im(E)<0). The approximate values for 
their positions obtained with the ACCC method and 4NF of equat ion (|13|) arc summarized in the last row of Table [I] 
The accuracy of these results is estimated to be of 50%. 
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FIG. 5: Comparison of resonance trajectories for ,V — 2~ tetraneutron, two different curves correspond calculations with 4NF 
eg. 1131 having length po=2.5 (dashed curve) and 2 fm (dot curve). The points correspond resonance positions for W being 
reduced from 36 in step of 4 MeV-fm -1 for po=2.5 fm curve and from 72 in step of 8 MeV-fm -1 for po=2 fm curve. 



In Table |U we have also displayed the strengths W' of the 4N force (|13[) at which the resonance trajectories crosses 
the imaginary-energy axis, slipping from fourth into third energy quadrant. One can argue that these values are pretty 
small and that a small correction of nuclear interaction (like the presence of attractive three-nucleon force) can push 
tetraneutron states back to the fourth energy quadrant (with positive real energy parts). On this point we would like 
to mention that the smallness of W' is only apparent and entirely due to the unrealistic long range character of the 
4NF we have chosen. W value would increase drastically if the range of the potential po is reduced to make 4NF 
more realistic . This fact is demonstrated in figure[5l where W value for 2 - tetraneutron state increases from 3.9 to 
10.8 MeV-fm -1 when po is reduced from 2.5 to 2 fm. This result shows that any realistic (p < 1.4 fm) multincutron 
force should be very strong to keep a multincutron resonances in the fourth energy quadrant. 

The J T = 2 + tetraneutron state represents an interesting case, since it shows the largest sensitivity to nn P-waves 
interaction. Some 3N and 4N scattering observables which are difficult to be reproduced with the existing models 
indicate a strong nn P-waves contribution. It has been suggested that these dis crep ancies in 3N and 4N scattering 
observables can be significantly improved by modifying nn P-waves within 20% |l5l llo . Il7l | . We have explored such 
a possibility and traced in fig. the tetraneutron resonance trajectory for Reid 93 interaction with nn P-waves 
enhanced by a factor 7=1.2. By doing so a slightly weaker critical strength of 40.02 MeV-fm -1 is required to bind 
tetraneutron, in comparison with 40.38 MeV-fm -1 of the original Reid 93 force. However, apart from a small shift 
in the E res (W) trajectory, such a modification of nn P-waves has not changed its qualitative behavior, ending up 
very close to its original value, always located in the third energy quadrant. These nn P-waves should be much more 
strongly enhanced, as much as creating dineutron resonances, to result in sizeable effects in tetraneutron resonance 
positions. 
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Finally, wc would like to remark that even if there was a resonance in the fourth energy quadrant having a small 
real energy part and a large imaginary one, it would be difficult to identify it experimentally. Resonance should have a 
rather small width T = —2Im(E res ) to produce a visible effect in the experimental cross section and a E = Re(E res ) 
centered Breit-Wigner shape. At most, it will give a weak enhancement in the cross section, hardly discernable from 
the background and not necessary centered around the E = Re(E res ). This makes very doubtful the perspective of 
physically observable tetraneutron resonances. Their eventual existence would imply a too strong modifications in 
the present nuclear Hamiltonians. 




FIG. 6: Sensibility of the 2 tetraneutron resonance trajectory with respect to nn P-waves. Solid line correspond to Reid 93 
nn interaction and dashed line was obtained with nn P-waves enhanced by a factor 7=1.2 



Our results are in qualitative agreement with the findings of Sofianos et al. [Hj , where authors were able to accurately 
determine the tetraneutron resonance positions in the third energy quadrant for positive parity states, although using 
S-wave MT I-III potential. Due to the small influence of P and higher nn partial waves on tetraneutron states, .S-wave 
models become very appropriate to study this system. 



IV. CONCLUSION 



Configuration space Faddccv-Yakubovsky equations have been solved with the aim of determining the positions of 
the four-neutron resonances in the complex energy plane. 

Realistic Reid 93 nn interaction model has been used. A systematic study of four-neutron resonances have been 
accomplished by first adding to the nuclear hamiltonian an attractive four-neutron force to artificially bind tetraneu- 
tron. The trajectory of the energy eigenvalue is then traced as a function of the strength of the additional force until 
it is fully removed. 

Two methods, namely Complex Scaling and Analytical Continuation in the Coupling Constant, were employed to 
follow these trajectories. 
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The low lying four-neutron resonance trajectories, corresponding to states with quantum numbers J 77 = , 1 , 2 , 
were shown to settle in the third energy quadrant (Re(E)< 0, Im(E)< 0)) well before the additional 4n force is 
completely removed. Furthermore, these resonances acquired a rather large imaginary energy T = 2Im(— E) w 20 
MeV and should hardly be experimentally observable. Tetraneutron compound - bound or resonant - can be created 
only in strong external fields and would disintegrate right after such a field is removed. 

Finally, we have demonstrated that the four-neutron physics is entirely determined by nn S- waves, namely -"-So one, 
which is controlled by the experimentally measurable nn-scattering length. All realistic nuclear interaction models 
should thus provide qualitatively identical results for tetraneutron resonances. This fact is supported by performing 
similar studies binding artificially tetraneutron states using various NN models |(|. 
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